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Abstract 

Data on the polarization observables T, P, and H for the reaction jp -A pir^ are reported. Compared to earlier data 
from other experiments, our data are more precise and extend the covered range in energy and angle substantially. The 
results were extracted from azimuthal asymmetries measured using a transversely polarized target and linearly polarized 
photons. The data were taken at the Bonn electron stretcher accelerator ELSA with the CBELSA/TAPS detector. 
Within the Bonn-Gatchina partial wave analysis, the new polarization data lead to a significant narrowing of the error 
band for the multipoles for neutral-pion photoproduction. 
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1. Introduction 

The measurement of the two double polarization ob¬ 
servables G [i| and E in photoproduction of neutral 
pions revealed significant differences between the data and 
the predictions from analyses such as MAID Q , SAID Q , 
and BnGa Q. Partly, large discrepancies were observed 
even at rather low photon energies. This was surprising 
since the reaction yp —s- pTr" is certainly the best studied 
photoproduction process. These discrepancies underline 
the importance of polarization observables for an interpre¬ 
tation of photoproduction data. 

In this letter, we report a measurement of further po¬ 
larization observables, called T, P, and H, for the reaction 

7P —>■ pTr" . (1) 

All three observables were determined simultaneously from 
the same measurement and provide the next important 


step toward a better understanding of tt" photoproduction. 
The target asymmetry T is a measure of the azimuthal 
asymmetry when the target nucleon carries polarization 
Pt in a direction perpendicular to the beam axis. P, often 
termed the recoil polarization observable, is a measure of 
the induced polarization of the recoiling nucleon. Here, 
P is determined from a double polarization measurement 
rather than from an experimentally more challenging di¬ 
rect measurement of the recoil polarization. This has the 
advantage that P can be determined in the very same mea¬ 
surement for almost the full solid angle, rather than by 
measuring P for specific points in angle and energy, as 
it has been done in the past. The observables P and H 
can be determined from azimuthal asymmetries using mea¬ 
surements with linearly polarized photons and transversely 
polarized target nucleons having polarization px. Part of 
the data presented here were used as a basis for an energy 
independent partial wave analysis (PWA) [^. They are 
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now included in the BnGa PWA and inultipoles for reac¬ 
tion 0 were determined. The multipoles were compared 
to those from MAID, SAID, Juelich2015 Q and earlier 
BnGa PWA-solutions. 


2. The experiment 


The experiment was performed at the Bonn Electron 
Stretcher Accelerator ELSA Q. Linearly polarized pho¬ 
tons were produced by scattering a 3.2 GeV electron beam 
off a diamond crystal [^. The crystal was oriented to po¬ 
sition the coherent edge at 950 MeV, leading to a polariza¬ 
tion maximum of = 65% at 850 MeV which declined to 
40% at 700 MeV. Two perpendicular settings of the beam 
polarization plane were used (named || and T). Photon 
ener^s were measured in a tagging system described in 
Ref. 0. 

The photon beam hits a butanol (G 4 HgOH) target with 
transversely polarized protons with a mean proton po¬ 
larization of pt ~ 75%. Data were taken with two opposite 
settings of the target polarization direction (named t and 

The incoming photons may produce a 7 r° in the re¬ 
action ©• The neutral pions were reconstructed from 
their 7 r° 2j decays in the Grystal Barrel (1320 GsI(Tl)- 
crystals) [Ilj and TAPS (216 BaFg crystals) [H, [l3| elec¬ 
tromagnetic calorimeters (see Fig. [T} which cover almost 
the full angular range down to 0 = 1 ° in the forward direc¬ 
tion. Protons from reaction o were detected in a three- 


layer cylindrical scintillation detector with 513 fibers [14 1 


surrounding the target, in 180 small organic scintillators 
in front of 90 forward GsI(Tl) crystals covering the an¬ 
gular range from 27.5° to 11.2°, and in organic scintil¬ 
lators mounted in front of each of the BaFg crystals. A 
CO 2 Cherenkov detector was installed in front of the BaFg 
crystals to identify background from electromagnetic reac¬ 
tions. The first-level trigger was derived from the tagger, 
the fiber detector, the forward calorimeters, and from the 
CO 2 Cherenkov detector as a veto; a second-level trigger 
used a FAst Cluster Encoder (FACE) [1^ and selected 
events with at least two distinct calorimeter hits in the 
full detector assembly. 

For further analysis only events with three distinct 
calorimeter hits were used. Adding events where only the 
two photons were measured would lead to an increased 
background contribution from other channels since sev¬ 
eral of the cuts discussed below can no longer be applied. 
For the events with three calorimeter hits, all three possi¬ 
ble combinations were treated as P 77 candidates, with the 
proton being treated as a missing particle. Kinematic cuts 
were applied to ensure momentum conservation. Examples 
for the missing mass and azimuthal angle difference distri¬ 
butions are shown in Fig. |31 Energy- and angle-dependent 
± 2 cr cuts were applied based on the corresponding distri¬ 
butions. In addition, a ±2(T-cut on the polar-angle differ¬ 
ence between the directions of the missing proton and the 
measured proton candidate was performed, and 7 p —?► p 7 r° 



Figure 1: (Color online) The central part of the detector system. The 
1320 CsI(Tl) crystals (blue and green) are read out via wavelength 
shifters and photodiodes (blue) or via photomultipliers (green), the 
216 BaF 2 crystals (yellow) in forward direction are read out with pho¬ 
tomultipliers. The target is surrounded by a 3-layer-scintillating fiber 
detector. To detect charged particles, scintillators are also placed in 
front of the forward CsI(Tl) (green) and TAPS crystals (yellow). 



Figure 2: The invariant 77 mass distribution for the butanol data, 
after all cuts discussed (see text) have been applied. The final data 
sample of events has a background contamination of less than 
1 %. 


events were selected by a zb 2 cr cut on the invariant 77 mass. 
Finally, a time coincidence was required between the tag¬ 
ger hit and the reaction products and a random-time back¬ 
ground subtraction was performed. This resulted in a final 
data sample containing a total of 1.4 million pir^ events. 
The background contamination was estimated from the in¬ 
variant 77 mass spectrum (see Fig. 13 ), assuming a linear 
behaviour of the background under the peak. It is less 
than 1% for all energies and angles. The selected events 
for each of the four combinations of beam and target po¬ 
larization directions were normalized to the corresponding 
photon flux and polarization degree for further analysis. 

Since a butanol target was used, not only polarized and 
unpolarized free protons contributed to the count rates but 
also reactions occurring on the bound unpolarized nucleons 
of the carbon and oxygen nuclei. Additional measurements 
using a carbon foam target were performed to determine 
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Figure 3: (Color online) (a,d) The missing mass distribution, with the proton as the missing particle, and (b,e) the azimuthal angle difference 
of TT^ and proton for reaction (1), for a invariant mass of VF = 1.46 — 1.48 (top) and 1.82-1.94GeV (bottom); butanol (□), scaled carbon 
(o), and the difference (A). The distributions are shown after all other cuts discussed in the text are applied. From these distributions, 
the dilution factor (c,f) is determined. The gray band indicates the systematic uncertainty in the dilution factor due to uncertainties in the 
determination of the carbon scaling factor s. Since only events with all three particles detected in the calorimeter are considered an acceptance 
hole occurs for large cos{0). The observed structures in the dilution factor are due to a combined effect of reduced efficiencies for clusters 
impinging onto detector boundaries and Fermi smearing. 


the so-called dilution factor d{W, cos 6) 


d{W, cos 9) 


N{ree{W,COs9) 

-^butanol(hf; COS O') 

-^butanol(hf; COS Vi3ound(hr; COS . 

A^butanol(Vr, COS 6>) 


which assumes that the nucleons bound in carbon and oxy¬ 
gen show the same response to the impinging photons. The 
carbon foam target had the same size as the butanol tar¬ 
get and approximately the same target area density as the 
carbon and oxygen part in the butanol. The carbon target 
replaced the butanol target in the frozen spin cryostat to 
match the experimental conditions of the butanol measure¬ 
ment as closely as possible. The flux-normalized carbon 
yield was compared to the flux-normalized butanol data 
using the distribution of the angle between the 7r° and pro¬ 
ton in the azimuthal plane outside the region where con¬ 
tributions from free protons can be expected (see Fig. [3]). 
Counting the yields for |A(() — 180°| > 20°, a global scaling 
factor s = 1.13 ± O.Olgtat. ± O.llgys. was determined based 
on which the dilution factor (Eq. [3]) was calculated using 


Vbound(bF, COs6>) = S • Vcarbon(bF,COs6»). (3) 


Taking the slightly different densities of the butanol and 
carbon targets into account, one would expect a scaling 


factor of s ss 1.1, which is in agreement with the value ob¬ 
tained from the data. Energy- and angle-dependent devia¬ 
tions from the global scaling factor s were investigated by 
determining s(lF, cos0) indepently for each bin and com¬ 
paring it to the global value. The observed deviations are 
of the same magnitude as their statistical uncertainty. The 
global value is used, with the deviation contributing to the 
systematic uncertainty. 

Figure |3] shows, for two energy bins, the azimuthal an¬ 
gle difference and the missing mass distributions of the 
butanol and scaled carbon data as well as their resulting 
difference. From these distributions the dilution factor as 
function of cos0 for each energy bin was determined. The 
dilution factor is quite large, around 0.9 at low energies and 
decreasing to around 0.6 at higher energies. Note that a 
dilution factor d = 1 corresponds to non-existent carbon 
background and therefore to no dilution. The observed 
structures (Fig. [3]) are due to a combined effect of reduced 
efficiencies for clusters impinging onto detector boundaries 
and Fermi smearing. They are reproduced in Monte Carlo 
simulations. At higher energies, the reduced missing mass 
resolution required a wider cut. Therefore the carbon con¬ 
tribution remaining after all cuts increased significantly, 
resulting in a smaller dilution factor. 
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3. The polarization observables 


In the coordinate frame of the detector system, with a 
being the azimuthal angle of the photon beam polarization 
plane in the || setting, jS the azimuthal angle of the target 
polarization vector in the t setting, and (j) the azimuthal 
angle of the produced 7 r°, the differential cross section can 
be written as 

^ (^) ’ ~ cos(2(q; - (t>)) + ptT sm{/3 - (j)) 

— P^PtP cos{2{a — 0))sin(/3 — (j)) 

+ P^PtH su\{2{a — (j)))cos{j3 — (j))}. (4) 


In a first step, the ordinary beam asymmetry for 
each bin in energy and angle was determined: 

1 jv I - iVii 

AiVbeam(<^) = — ' ' cos( 2 (a - 0 )). (5) 


P'1 N± + IV|| 


A typical distribution is shown in Fig. 0^. The resulting 
beam asymmetries agree very well with previously re¬ 
ported measurements [Ol, [la, ll7| although the data sample 
contains in part reactions off nucleons of the C/O-nuclei. 
Results for are not shown here. 

To determine T, the asymmetry of the data sets with 
respect to the target polarization directions t and J, was 
used, resulting in 


A^target(</>) = ^ = T • sin(/3 - <^) ( 6 ) 

with d given by Eq. 0] The target asymmetry was deter¬ 
mined by a fit to the AA^target distributions as shown 
by the example in Fig.|3)D. Figure [5] shows AA^target(<^) for 
three different angular bins (VF=I.82-I.94GeV). These 
distributions underline the strong dependence of T on the 
scattering angle 9. 

P and H can be extracted from the data by considering 
the linear beam-polarization plane (|| and _L) in addition 
to the target polarization, leading to 


AA^beam-target 


1 - Nu) - (A^llt - ^llt) 

dp'/PT (A_l|- 1-Af_L4_)-I-(A^||.|.-I-Af||j_) 
Psin(/3 — (p) cos(2(a — (p)) 

+ H cos{f3 — (p) sin{2{a — <p)) (7) 


with average beam polarization pj and a = 45°, the di¬ 
rection of the polarization plane in the || setting. Again, 
the observables could easily be determined by a fit to the 
AA^beam-target((/') distributions, see Fig. 01:. 

The AN{(p) distributions were fitted for each bin in 
energy W and angle 9. The resulting confidence level (CL) 
distribution for the fits is flat with a distinct increase well 
below 0.1%, as can be seen in Fig. 01 This increase was 
traced to fits of distributions with very low statistics (due 
to a low cross section or low acceptance). These fits were 
very sensitive to background fluctuations. These data were 


excluded from further analysis by performing a CL-cut at 

0 . 1 %. 

The data on T, P, and H as functions of center-of-mass 
energy W and angle cos0 are shown in Figs. [7]-[ni All ob¬ 
servables exhibit a strong angular dependence which also 
changes significantly with W. The systematic uncertainty 
shown includes the uncertainty in the degree of photon 
(4%) and proton (2%) polarizations, in the dilution factor 
(l%-4%), and an additional absolute uncertainty (0.01) 
due to the remaining background contribution. For fur¬ 
ther details on the estimation of the systematic uncertain¬ 
ties see Refs. 

The polarization observable P describes the polariza¬ 
tion of the outgoing proton in the direction perpendicu¬ 
lar to the scattering plane. Here, it was determined indi¬ 
rectly from the correlation of beam and target polariza¬ 
tion. Thus, both observables, P and H, were measured 
only in the energy region in which the photon beam car¬ 
ried a significant linear polarization; the results are hence 
restricted to the 665 to 930 MeV photon energy range. The 
target asymmetry T, shown in Fig. [71 does not require po¬ 
larized photons, and the data allowed us to determine T 
up to = 1900 MeV. Above this energy, the count rates 
were small, and we do not present those results here. 

The new dat a ag ree well with previously reported mea¬ 
surements of P [l9| and T but exceed the old data in 
precision and coverage in angles and energy. For H no ear¬ 
lier data exist in this energy range, the older data [2]J are 
limited to > 1300 MeV. Our data up to 930 MeV were 
used as a basis for an energy independent PWA Q. The 
high energy T-data are presented here for the first time. 
The new data sets have been included in the BnGa-PWA 
as discussed in the following. 


4. Partial wave analysis 


The data were fitted within the BnGa multi-channel 
partial wave analysis. Gompared to our solutions BnGa2011- 
01 and BnGa2011-02, further data were included [l|,|2, 25|. 
Relevant for the yp —> p 7 r° multipoles are the new data 
on T, P, and H as well as our recently published data 
on G [l| and on P 0. Figures [TUI and [111 compare the 
newly determined multipoles with those of BnGa2011-01 
and BnGa2011-02. The error bands for the BnGa2011 so¬ 
lutions were derived from the (Icr) spread of 12 different 
solutions with different assumptions on the ingredients: 
the number of poles in the = 2>/2'^ wave was 3 or 4, 
in the = 5/2+ wave 2 or 3, the Af(1700)3/2“ width 
of the pole converged to a wide (~600MeV) or a narrow 
(^^250 MeV) value, a K-matrix formalism was used or, al¬ 
ternatively, an A^/P-parametrization [^. The ranged 
from its minimum to y^j^ + 800. (Note that the ab¬ 
solute y^ value is meaningless since part of the data are 
multiparticle hnal states and fitted in an event-based like¬ 
lihood fit. The log likelihood value is then converted into 
a pseudo-y^ Q.) 
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Figure 4: (Color online) Yield asymmetries as a function of (p for the energy bin 1.46 GeV< W < 1.48 GeV. (a) AYbeam(‘?^)! (t>) A Ytarget (<?!>), 
(c) AYbeam-target('?^) fitted by the function of Eq. Eq. and Eq. Q, respectively. 



(|) 0 (j) 


Figure 5: (Color online) Yield asymmetry AY^target as a function of 0 for three angular bins in 9 in the 'yp invariant mass window 1.82 GeV 
< VF < 1.94 GeV. 


In the new fits, we started from the same solutions 
and re-optimized all parameters. All fits converged, but 
6 fits resulted in a larger than the new 
(Xnewmin) 1000 units Or more. These fits, mostly those 
with only 3 poles in the = 3/2''" wave, were then re¬ 
moved from the error analysis. The resulting error bands 
for all remaining solutions within Xnew min + 
shown in Figs. [TUI and [IT] The x^-range for the new solu¬ 
tions was chosen consistently with the BnGa2011 solutions 
for which the multipoles have first been shown in 241. 

The new error bands are significantly smaller than the 
previous ones. Averaged over all multipoles and ener¬ 
gies, the errors are reduced by a factor of 2.25. Exam¬ 
ples for multipoles which are substantially better defined 
are the multipole to which the Roper resonance and 
V(1710)l/2“ contribute and the multipole to which 
the = 5/2“-resonances contribute. In most cases, the 
old error bars cover the range of new solutions, the solu¬ 
tions are at least compatible with each other at the 2 ct 
level over the full mass range. 

There are few changes to the multipoles only: The Ef 
multipole leading to the excitation of N and A resonances 
with = 3/2+ has kept its structure but its imaginary 
part has increased in strength in the fourth resonance re¬ 
gion. Similar changes are also visible in the real part. The 
real part of the multipole (J^ = 5/2“) has become 
smaller in absolute value in the W = 1500 MeV region. 


while in the imaginary part changes are observed in the 
high mass region around W = 2100 MeV. In all cases, 
where discrepancies with the MAID and SAID solutions 
were observed, these discrepancies remain while the con¬ 
sistency between BnGa2011 and BnGa2014 is rather good. 
The multipoles from Juelich2015 Q show also significant 
differences compared to BnGa2014. 

Summarizing, we have reported a simultaneous mea¬ 
surement of the polarization observables T, P, and H. 
With the data presented here an additional step toward a 
complete experiment in 7r®-photoproduction off the proton 
has been made. The data provide a more precise determi¬ 
nation of the photoproduction multipoles governing the 
photoproduction of single neutral pions off protons. 
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Figure 7: (Color online) The polarization observable T as function of the invariant mass W (in GeV) and of the scattering angle cos0. 
The systematic uncertainty is shown as dark gray band. An additional systematic uncertainty on the photon energy (from =6.5 MeV at 
the lowest to 2.3 MeV at the highest energy bin) is not shown. The low energy data were presented in [^. Earlier data (gray, (red, online)) 
are from [Soll . The solid line represents our best fit BnGa2014. The data are compared to predicitions (dashed curves) from BnGa2011-02 
(black), MAID (light gray, (green, online)), SAID CM12 (dark gray, (blue, online)), and Juelich2015 Q (gray, (magenta, online)). 
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Figure 8: (Color online) The polarization observable P as functions of the 7 p invariant mass W (in GeV) and of the scattering angle cos0. 
The systematic uncertainty is shown as dark gray band. An additional systematic uncertainty on the photon energy (from =6.5 MeV at 
the lowest to 5.4 MeV at the highest energy bin) is not shown. Earlier data (gray, (red, online)) are from [T^ . The solid line represents our 
best fit BnGa2014. The data are compared to predictions (dashed curves) from BnGa2011-02 (black), MAID Q (light gray, (green, online)), 
SAID CM12 (dark gray, (blue, online)), and Juelich2015 a (gray, (magenta, online)). 



Figure 9: (Color online) The polarization observable if as functions of the 'yp invariant mass W (in GeV) and of the scattering angle cos0. 
The systematic uncertainty is shown as dark gray band. An additional systematic uncertainty on the photon energy (from =6.5 MeV at 
the lowest to 5.4 MeV at the highest energy bin) is not shown. The solid line represents our best fit BnGa2014. The data are compared to 
predicitions (dashed curves) from BnGa2011-02 (black), MAID (light gray, (green, online)), SAID CM12 (dark gray, (blue, online)), 
and Juelich2015 0] (gray, (magenta, online)). 
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Figure 10: (Color online) Multipole decomposition of the 'yp —>■ transition amplitudes, imaginary part. The light (red) shaded areas give 

the range from a variety of different fits derived from solution BnGa2011-01 and BnGa2011-02 |^. The dark (blue) shaded area represents 
the range of solutions when the new data are included in the fit. The black curves represent the MAID fit [^, the light (green) solid curves 
SAID-CM12 [^, the light (green) dashed curves SAID-SNll i , and the magenta curve the Juelich2015 Q solution. For the BnGa-multipoles 
an error band (see text) has been determined. Such an error band is presently not provided by the other analyses. 
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Figure 11: (Color online) Multipole decomposition of the 7 p —>■ transition amplitudes, real part. The light (red) shaded areas give the 
range from a variety of different fits derived from solution BnGa2011-01 and BnGa2011-02 Q. The dark (blue) shaded area represents the 
range of solutions when the new data are included in the fit. The black curves represent the MAID fit [^, the light (green) solid curves 
SAID-CM12 [^, the light (green) dashed curves SAID-SNll i , and the magenta curve the Juelich2015 0 solution. For the BnGa-multipoles 
an error band (see text) has been determined. Such an error band is presently not provided by other analyses. 
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